Numerical heterodyne holography with two-dimensional photodetector arrays 
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We present an original heterodyne holography method for digital holography 
that relies on two-dimensional heterodyne detection to record the phase and 
the amplitude of a field. The technique has been tested on objects as much 
as 13 mm in size. Consistency checks were performed, and high-resolution 
images were computed. We show the requirement for a spatial filter to select 
properly sampled near-axis photons. Heterodyne holography is superior to 
off-axis digital holography for both field of view and resolution. 
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As was demonstrated by Gabor^, in the early 1950. s, 
the purpose of holography is to record the phase and 
the amplitude of the light coming from an object un- 
der coherent illumination. Classical holography does 
not provide straightforward access to the holographic 
data. For quantitative analyses of those data, in digi- 
tal holograph}*^ (DH) photographic films were replaced 
by two-dimensional electronic detectors. In both digital 
and thin-film holograms a ghost field and the remaining 
part of the reference field are superimposed upon the re- 
constructed object field^. A solution to this problem is 
to tilt the reference beam with respect to the object^, 
to separate physically the spatial-frequency components. 
Although this off-axis technique is acceptable for high- 
resolution holographic films, it is hardly compatible with 
the limited resolution of digital holograms. All these ar- 
tifacts arise from the fact that among all these meth- 
ods only one field quadrature is measured. Measuring 
both quadratures requires recording at least two inter- 
ferograms with distinct reference phases^. In classical 
holography this is done with a thick plate that records 
several intensity fringes in depth. 

In this Letter we describe a heterodyne holographic 
scheme in which the reference beam is dynamically phase 
shifted with respect to the signal field. This shift pro- 
duces time-varying interferograms on a two-dimensional 
sensor. In our experiment the phase shift is linear in 
time (frequency shift). Intensity / in the detector plane 
results from the interference of the signal field with the 
(5/-shifted reference field: 



I{t) = \Es + ERexp{2t7Tdft)f 



(1) 



the signal and the reference fields, respectively. L in- 
tensity /; {I — 0...L — 1) measurements are performed 
within a Sf period at t; — 2-111/6 f. We obtain Es by 
demodulating / : 
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where * is the complex conjugate. For L = 4, Es is pro- 
portional to {Iq— l2) + i{h ~h)- Heterodyne holography 
(HH) thus measures the phase, using the information ob- 
tained at different times, and DH extracts the phase from 
measurements made of different pixels^. In both cases the 
sampling theorem^ restricts the largest admissible field- 
of-view angle to 



1^1 (^max — ^/ {'2dpixel) 



(3) 



where A is the wavelength and dpi^ei is the pixel spacing. 
In DH the object must be off axis, which yields the con- 
straint [3] that 9 > 9min = N dpixei / {2D) , where N is the 
number of pixels and D is the sensor-to-object distance. 



The DH field of view, 9,, 
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is thus much smaller 



where £^5 and En represent the complex amplitudes of 



than for HH: 16*1 < 9max- Whereas 9min < 9max- restricts 
DH to the far field with respect to Dmin = Ndpixei /A, HH 
works in both the near field and the far field modes. In 
the far field, both HH and DH angular resolution reaches 
the diffraction limit, 9diff = X/{Ndpixei), which corre- 
sponds to the two dimensional sensor's size. In HH the 
number of resolved pixels (i.e., field of view/angular res- 
olution) is N. This remains true in the near-field regime. 

Although HH works on ref lection, we focus here on 
a transmission configuration. The general setup, shown 
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Fig. 1. Experimental setup: L, He.Ne laser; other abbre- 
viations defined in text. 
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backtransforms the field to real space. The SF, 02, the 
BS, and the CCD are kept in alignment, so the selected 
photons reach the CCD nearly parallel to the reference 
beam. Equation [3] yields a rectangular SF of dimen- 
sions dr^^y = 2taii(X/2dpixei^ y)f2 {dx = 1-84 mm and 
dy — 1.90 mm). The 01-02 optical system magnifies the 
incoming beam by a factor /1//2 = 1/2. Our setup thus 
has X and y fields of view of ±1.05° and ±1.09° , cor- 
responding to an equivalent CCD with magnified pixels 
d'pixei = '^dpixei- 01 is fixed, and the SF, 02, the BS, 
and the CCD camera can be x and y translated by step 
motors. When the SF is not centered on the optical axis 
of 01, the system records holograms that correspond to 
tilted k components, as depicted in Fig. 2. Those dis- 
placements allow us to center the SF accurately. This 
feature will be used in the future to record wider-field-of- 
view holograms by merging several holograms with dif- 
ferent fc-component tilts. 

To test our setup quantitatively, we recorded the holo- 
gram of a narrow slit 9.5 mm high and w = 35/im wide 
at = 60 cm with a 480-ms integration time. We 
got Es{x,y,z = 0) on the equivalent CCD. To compute 
Es{x,y,z = D) on the slit we considered the Fresnel 
propagation of Es from z = to z = D, which can be 
formally expressed as an x and y convolution product 
(symbol 0): 



Esix,y,D) = e'^°''[Piz = D)(^^Esix,y,0)] 



P{x,y,z) = TT-exp 
iXz 



(4) 



Fig. 2. SF system: LO, local oscillator; other abbrevia- 
tions defined in text. 



in Fig. [TJ is a Mach.Zehnder interferometer composed of 
two beams (reference and signal) from the same coher- 
ent source (a 5-mW He.Ne laser). The two beams are 
expanded by beam expanders BEl and BE2 and com- 
bined by a beam splitter (BS) on a CCD camera with 
Nx = 768 and Ny = 576 pixels, dpixei^ — 8.6 fim and 
dpixeiy = 8.3 Ilia. All X and y subscripts below refer to x 
or y axes, respectively. We 6f shift the reference beam 
by combining two acousto-optic modulators, AOMl and 
A0M2, working at Af + 5f and —A/, respectively, with 
A/ = 80 MHz. 6f = 6.25 Hz is equal to one quarter 
of the CCD image frequency // = 25 Hz {L = 4). The 
video frame is acquired by an 8-bit analog frame grabber 
(Matrox Meteor), and a Pentium II 450-MHz computer 
calculates the complex field in real time. 

To select the signal on- axis photons, we insert between 
the object and the BS a spatial filter system (Fig. 2) 
composed of two confocal objectives, 01 (focal length, 
/i = 50 mm) and 02 (/2 = 25 mm), with a spatial fil- 
ter (SF) in their common focal plane. 01 transforms 
the field into its fc-space components in its focal plane 
where the SF selects the photons that fulfill Eql31 02 



We obtained a 13.2 mm x 9.6 mm (equivalent CCD 
size) contact image of the slit in its plane [Fig. EJa)]. 
The z = field is calculated on a 1024 x 1024 grid 
by bilinear interpolation of the 768 x 576 experimental 
points, and the convolution product is calculated by the 
fast-Fourier-transformmethod^, with the grid size kept 
constant over z. Figure [HJb) shows the intensity along 
the A A' cut. Points are experimental data on the fast- 
Fourier-transform grid. The solid curve is the calculated 
diffracted intensity of an ideal wide slit (dashed lines). 
Both diffraction and spatial averaging over the magni- 
fied pixels are calculated; the agreement is good. D is 
larger than Dmm^ = 35.9 cm and Dmin = 25.0 cm, so 
diffraction is the main contributor to resolution; here it is 
assumed to be equal to the distance that yields a sharp 
contact imaging. Comparing D with the 01-to-object 
distance, we find that the equivalent CCD is located ^ 27 
cm behind 01. This location does not depend on D and 
agrees with the positions of 01, 02, and the CCD. 

Figure [3] shows the contact image of a transmission 
U.S. Air Force test chart target at D = 58 cm. At 
the right. Fig. 2] is exactly on focus, whereas at the 
left the target is 1.2 mm out of focus. The blur is 
due to a defocus equal to the x depth of focus (DOF): 



DOFx,y = A(l 



NAiy^^ 



f/NAlyf [DOFx = 1.22 mm, 



DOFx = 2.32 mm). Here NAx,y = Nx,yd[^,,,i^ J D is 
the numerical aperture. 
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Fig. 3. (a) Contact image of the slit (b) cut along AA': 
filled squares, experimental intensity points; dashed lines, 
initial slit; solid curve, theoretical intensity profile. 



Figure \5\ shows the central part of the contact images 
of black characters printed on a diffusive plastic sheet at 
D = 30 cm with (top) and without (bottom) a SF. Here 
the plastic substrate scatters light to angles larger than 
Smax ■ Without a SF the off-axis photons scramble the 
useful information, reducing the image contrast. This 
phenomenon is a spatial- frequency- folding effect. 

As was shown in this Letter, HH measures the field 
within the SF sampling cone without information loss. 
HH has many other advantages. By measuring the field 
in a z = plane we are able to compute it at any z plane, 
thus getting three-dimensional information. Because the 
detected interference arises from the heterodyne mixing 
of a strong reference with a weak signal, HH is expected 
to be sensitive up to the photon level. Compared with 
intensity imaging, HH offers a huge dynamic range be- 




Fig. 4. Contact image of a U.S. Air Force test chart 
target a.t D — 58 cm. Right, image on focus; left, 1.2- 
mm out-of-focus image. 



cause it detects amplitude. The time-varying character 
of HH provides good immunity against any dc and non- 
harmonic (with respect to Sf) errors. The CCD image 
phase shift (90°) is highly accurate because it results from 
a frequency-offset Sf. It is also possible to shift the het- 
erodyne frequency Sf for sideband heterodyne detection. 
Making Sf = (1/ 4:)fi+fmod allows the vibration of an ob- 
ject at fmod to be detected. HH works with a low-power 
cw laser. HH may be performed with low coherent laser 
sources to select a narrow space slice on the object where 
coherence is conserved. 

Within the numerous possible applications of the 
method, we have explored diffusing medial and aper- 
ture synthesis applications, and further experiments are 
in progress. In both cases, a SF is useful. 
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Fig. 5. Contact images of black characters printed on 
diffusive plastic sheets with (top) and without (bottom) 
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